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I describe a recently derived stochastic approach to infla- 
ton dynamics which can address some serious problems asso- 
ciated with conventional inflationary theory. Using this the- 
ory 1 derive an exact solution to the stochastic dynamics for 
the case of a \(f) 4 potential and use it to study the generated 
primordial density fluctuations, ft is found that on both sub 
and super-horizon scales the theory predicts gaussian fluctua- 
tions to a very high accuracy along with a near scale-invariant 
spectrum. Of most interest is that the amplitude constraint 
is found to be satisfied for A ~ 10 -5 rather than A ~ 10 -14 of 
the conventional theory. This represents a dramatic easing of 
the fine-tuning constraints, a feature likely to generalize to a 
wide range of potentials. 



The inflationary universe scenario asserts that, at some 
very early time, the universe went through a de Sitter 
phase expansion with scale factor a(t) growing as e Ht . 
Inflation is needed because it solves the horizon, flatness 
and monopole problems of the very early universe and 
also provides a mechanism for the creation of primordial 
density fluctuations. For these reasons it is an integral 
part of the standard cosmological model [1] . 

The inflationary phase is driven by a quantum scalar 
field with a potential V($>), that can take on many dif- 
ferent forms that satisfy the 'slow roll' conditions. In the 
conventional approach to inflaton dynamics [1], the in- 
flaton field $ is first split into a spatially homogeneous 
piece and an inhomogeneous piece 



(tf) = (27T) 



' 2 Hh. 



(3) 



$(x,s) = (j)(s) + ip{x,s). 



(1) 



The dynamics of the <j> is then postulated to obey the 
classical 'slow roll' equation of motion 



3H 



= 0. 



(2) 



This equation governs the dynamics of <fi which drives the 
inflationary phase. It is also possible to discuss the gener- 
ation of primordial density fluctuations using ip. Assum- 
ing that p 3> ip, it can be shown that ip is described by 
a free massless minimally coupled quantum scalar field. 
During exponential inflation the quantum fluctuations of 
ip grow as [2] 
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These quantum fluctuations are then identified with the 
classical fluctuations which generate primordial density 
fluctuations [1,3]. It is important to note here that inter- 
actions between the coarse-grained field p> and its fluctu- 
ations ip are ignored. This is possible in this approach 
because the density fluctuations are directly identified 
with (ip 2 ). 

Consistent with the conventional approach above is the 
'Stochastic Inflation' program initiated by Starobinsky 
[4] and further developed by others [5]. In this case the 
field <p obeys 

• V'U) # 3/2 , s 

where F w (t) is a zero mean gaussian white noise source 
of unit amplitude. In this case p describes the field $ 
coarse-grained over a volume determined by the de Sit- 
ter Hubble radius. We will refer to <p as the local order 
parameter. In this method the observable universe is 
comprised of many patches each with its own local order 
parameter whose dynamics obeys (4). Spatial inhomo- 
geneities arise because the local order parameter in each 
patch can take on different values by virtue of the noise 
in (4). Equation (4) has been the basis for many appli- 
cations including studies of the generation of primordial 
density fluctuations [6,7] and the very large scale struc- 
ture of the universe [8] . 

A problem with the conventional approach is that it is 
assumed, without justification, that the local order pa- 
rameter <p can be treated as a classical order parameter, 
and that the quantum fluctuations of ip are equivalent to 
classical fluctuations. Since cj> and ip are treated as de- 
coupled closed quantum systems it is impossible for this 
method to explain the quantum-to-classical transition of 
(p and ip from first principles. Another more serious prob- 
lem comes from directly identifying the quantum fluctu- 
ations (ip 2 ), with the classical fluctuations that generate 
primordial density fluctuations. This scheme leads to an 
overproduction of primordial density fluctuations which 
can only be avoided by unnaturally fine-tuning the cou- 
pling constants in the inflaton potential. This is the well 
known fine-tuning problem of inflation. 

Several authors have previously suggested that these 
problems arise because the conventional approach to cal- 
culating primordial density fluctuations is inconsistent 
with the established methods of non-equilibrium statisti- 
cal physics. This was first pointed out by Hu and Zhang 
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[9] and further developed in [10] (see also Lombardo and 
Mazzitelli [11] and Morikawa [12]). Morikawa [13] first 
suggested that this inconsistency was the origin of the 
fine-tuning problem. Calzetta and Hu [14] and more re- 
cently Calzetta and Gonorazky [15] independently and in 
much greater detail addressed this issue for a A</> 4 theory. 

While the conventional approach may be the only pos- 
sible one for a free field, in [16] an alternative has been 
developed for interacting fields which does address the 
problems outlined above. The theory is similar in style to 
the conventional stochastic inflation program but differs 
in a fundamental way. In this theory we no longer iden- 
tify the quantum fluctuations (V> 2 ) directly with the clas- 
sical fluctuations that generate primordial density fluctu- 
ations. The new role of the field ip is to provide a noise 
source (via backreaction) in the quantum dynamics of 
the local order parameter <fi. This is nothing but an ap- 
plication of the well known quantum Brownian motion 
paradigm of non-equilibrium statistical physics (see [17] 
and references therein). The field ip plays the role of 
an environment which couples to the system <j) and in- 
directly generates fluctuations 5(f) in the system. This 
environmental noise will generate quantum decoherence 
which is the process that leads to entropy generation and 
the quantum-to-classical transition of the order parame- 
ter and its fluctuations. We then identify the resulting 
classical fluctuations of the local order parameter 5(j), as 
those which lead to density fluctuations, rather than the 
quantum fluctuations derived directly from (f/> 2 )- Clearly 
in this approach the interaction between 4> and ip is criti- 
cal. As well as addressing the quantum-to-classical tran- 
sition problem, this theory leads to a dramatic easing of 
the fine tuning constraints, a problem that has plagued 
the conventional approach to inflaton dynamics. 

In the classical limit of this theory, the dynamics of the 
local order parameter is described by [16] 



4> + SH4> + v'{4) = §^v'"{<j>)F c (t) 



(5) 



where F c is a colored gaussian noise of unit amplitude 
with a correlation time of the order H . The origin 
of the noise is the backreaction of quantum fluctuations 
with wavelengths shorter than the coarse-graining scale. 
The noise correlation function is ultraviolet finite and 
also turns out to be independent of any ultraviolet cut- 
off. The noise is of a multiplicative nature because its 
origin is the mode-mode coupling induced by the self- 
interaction of the inflaton. For a free field the stochastic 
term vanishes. This is because the environment if) and 
the system <j> now decouple and the conventional situ- 
ation is recovered. Significant simplification of (5) was 
obtained by invoking the standard slow roll assumptions. 
This made it possible to show that neglecting the poten- 
tial renormalization and non-local dissipation terms was 
a good first order approximation in the early slow roll 
phase. 



The coarse-graining scale must be greater than the 
Hubble radius. This condition allows us to ignore the 
spatial gradient term in the system sector. The theory is 
essentially independent of the coarse-graining scale when 
this condition is met. This robustness to the nature of 
the coarse-graining is an important virtue of the theory. 
We also ignore information about spatial correlations be- 
tween the order parameters of different regions. This al- 
lows a description based on a single degree of freedom. 
Using the slow roll assumptions we can drop the incrtial 
term in (5) and approximate the colored noise by a white 
noise. In this case we find that (5) becomes 



3H 



#1/2 



■V'"^)F w (t). 



(6) 



F w is a white noise of unit amplitude that is interpreted 
in the Stratonovich sense (since it is an approximation to 
a colored noise). We also interpret the noise in (4) the 
same way, though in this case one is also free to use the 
Ito interpretation. Equation (6) is the result we will use 
in this letter. 

Equation (6) is valid only as long as the slow-roll ap- 
proximation is valid. This approximation is valid when 
the slow-roll parameters e and \r)\ [18], which are defined 

by 



and 



<4>) 



rfa ( V\4>) 
16tt ^ V(</>) 



v"{4>) i (YWX 
vtf) 2 V v(<f>) J 



(7) 



(8) 



are both less than 1. During inflation we have e(<j>) < 1, 
and the local order parameter rolls down the potential 
hill according to (2) from its initial value (f>o- Inflation 
ends at a field value <f) e which is determined by 



e(0e) 



V 



16tt V V(cp e ) 







= 1. 



(9) 



At this point slow-rolling ends and the reheating phase 
commences. Our discussion in this letter will be re- 
stricted to the potential 



V{cj)) = Xcj) 4 /4 



(10) 



for which equation (6) will be exactly solvable. For the 
potential (10) we find from (9) that (j) e = 0.56ra p /. This 
sets a lower bound on </>. The second slow roll condition 
is always satisfied up to this value. 

The number of e-folds of inflation which occur when 
the field evolves from <p to <p e is [1] 
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Smoothness on scales comparable to the current observ- 
able universe requires N > 60. This places a lower limit 
on the initial field value </> > 060, where N((f> 60 ) = 60. 
For the potential (10) we find from (11) that 



N{<l>) = - T {4> 2 -<j>l). 



in 



(12) 



pi 



From this we can deduce 4>§o = 4.4m p /. In this letter, as 
is common in most inflation models, we assume that the 
observable universe leaves the horizon during inflation 
at 60 Hubble times before the end of inflation, i.e when 
the inflaton field has the value 4>qo- The smallest scale 
for which density fluctuations can be probed (~ IMpc) 
will leave the horizon about 10 Hubble times after the 
observable universe. From (12) we find this corresponds 
to the field value </> 50 = 4.0m p ;. Significant is that the 
critical field values <f> e and <^60-50 are independent of A. 
The initial field <f>$ cannot be arbitrarily large. It must 



where (b p i is the field at the Planck 



V 



While 6 P and 



also satisfy <po < 
boundary which satisfies V(<b p i) 

050-60 are independent of A, this is clearly not so for (b p i. 

In this letter we wish to calculate the statistical proper- 
ties of observable density fluctuations predicted by (6) for 
the potential (10). Our aim is to demonstrate that this 
new approach to inflaton dynamics is consistent with the 
observed near gaussian and scale invariant density fluc- 
tuations. It is well known from gauge-invariant analysis 
that the amplitude of a density fluctuation that crosses 
back inside the horizon after inflation, can be deduced 
from the quantity [1,3] 



Sp 



,H{(b) 



(13) 



evaluated at the time the fluctuation scale of interest 
crossed outside the horizon during inflation. The power 
spectrum A(fc) is related to the mean square density fluc- 
tuations Sp/ p via 



Sp 
P 

from which we obtain 



A(k) 



A(k)d\nk 



Sp 



d\nk \ p 



(14) 



(15) 



The power spectrum is just the contribution, at a given 
time, to the mean square fluctuations generated in a Hub- 
ble time. We are interested in the power spectrum over 
observable scales. Therefore we will evaluate the RHS of 
(15) at 4>qq which is the classical field value at the time 
the scale of the observable universe leaves the horizon 
during inflation. It is usual to assume that the density 
fluctuation power spectrum A(k) can, within the observ- 
able range of k, be written as 



A(Jfe) = Ak 



n-l 



(16) 



where A and n are the amplitude and spectral index of 
the density fluctuations. For n = 1 we have a scale- 
invariant power spectrum of density fluctuations. From 
(16) we find that the spectral index can be determined 
from the power spectrum by 



n = 1 + 



1 dA(k) 
A(jfc) dhik' 



(17) 



Clearly our first task is to calculate (13) for which we 
need to solve (6) in order to obtain 5(b. As is commonly 
done in studies of the Starobinsky equation (4) [6,7], we 
include the effects of backreaction by simply assuming 
that H is slowly varying as 



H\(b) 



8tt 



V{(b). 



(18) 



This is possible since H changes little over a Hubble time 
clue to the slow rolling of the inflaton. Upon substituting 
(10) and (18) into (6), and changing to the dimcnsionless 
variables x — (b/m p i and r = tm p i, we obtain 



dx = —fxdr + gx 2 F w {r)dT 



where 



(19) 



(20) 



With the new variable x = — 1/a; we find that (19) be- 
comes 




2tt 



1/4 



— A 5 /4. 
247T 3 V 3 



dx = fx + 9F w (r). 



(21) 



This equation now describes an Ornstein-Uhlenbeck pro- 
cess for which the solution is 



X(r) 



= Xoe fT +ge fT f e-f s F w (s)ds. (22) 
Jo 



Therefore our solution to (19) is 

x c (t) 



x(t) = 



(23) 



1 -gx J T e f s F w (s)ds' 

where the classical deterministic solution x c (t) is 

x c (T)=x e^ T (24) 

with xo as the initial field value. To obtain the fluctua- 
tions 5(b we make a gaussian approximation 

x(t) ~ x c (t) (l + gx J e-f s F w (s)d S ^j (25) 

to (23) which will be justified later. From this we obtain 
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{5xf = {{x{r)-x c {r)f) 
A 2 



24tt 5 



xl{r)(xl-xl(r)). 



(26) 



We can compare this result directly to that predicted 
by the Starobinsky equation (4). Upon substituting (10) 
and (18) into (4), we find the exact solution to the 
Starobinsky equation is 



x(t) 



x c {t) 



l-2fcr§/ T e- 2 f s F w (s)ds 

x c (t) 

1 -hx\ll e-^ s F w {s)ds ' 



where 



1 /2ttA\ 
h= ^{ — ) 



3/4 



(27) 



(28) 



Making the same gaussian approximation as previously 
we find that 



{5x s f^ A x 2 (t)(x 4_ x 4 (t)) 



(29) 



where we use the s subscript to denote quantities evalu- 
ated from the solution to the Starobinsky equation. 
Using (26), (18) and (24) we find that (13) becomes 

(7) 2 = i A2x ' (r) ^" x ' (r) )- (30) 

To calculate the spectrum (15) we need to relate the clas- 
sical field value at some time to the scale k that is crossing 
the horizon at the same time. We do this by first con- 
sidering the number of e- folds N(k) between the horizon 
crossing of a scale k and the end of inflation. We are as- 
suming that the scale of the observable universe left the 
horizon 60 Hubble times before the end of inflation. We 
can therefore write N{k) as [1] 



N(k) = 60 + ln(fc*/fc) 



(31) 



where fc* is the scale of the current observable universe. 
From (31) and (12) we find that 



x 2 (r) = - [60 + ln(fc./fc)l 

7T 



xl 



(32) 



This determines the value of the field at the time a scale 
k crosses the horizon. This is the result we will use to 
calculate the spectrum. Substituting (30) into (15) and 
using (32) we find 



A(k) = -^x 2 c (T)(3x 2 c (r)-2x 



(33) 



The explicit k dependence of this function is contained in 
x c (t) via equation (32). The spectrum has only a very 



weak logarithmic k dependence. This will clearly give 
rise to a near scale-invariant spectrum. We can calculate 
the spectral index from (33) by using (17) and setting x c 
and xq to xqq. Observable constraints on the amplitude 
of density fluctuations require setting the LHS of (33) to 
10~ 10 . Setting x c and x to x 6 o allows us to calculate a 
value for A. The results are 



1.2 x 10" 



n - 



0.93. 



(34) 



We can perform an exactly analogous calculation using 
the fluctuations (29) derived from the exact solution of 
the Starobinsky equation. We find the familiar result 



A, = 6.6 x 10~ 15 , 



0.92. 



(35) 



We see that equation (6) leads to a dramatic easing of the 
fine-tuning required, but still gives a near scale-invariant 
spectrum of density fluctuations. Similar easings of the 
fine-tuning constraints have been reported in [13-16]. 

Also of great interest is a measure of how much the 
probability distribution of fluctuations deviates from 
gaussian. When the magnitude of the stochastic term 
in the denominator of (23) is small (<C 1), it becomes 
possible to make the gaussian approximation (25) to the 
exact solution (23). We are therefore led to define 



d = 



X X 



(36) 



as a simple measure of the deviation from gaussian. (The 
denominator \/2f of (36) is the typical fluctuation size 
of the stochastic term in (23) in the stationary limit.) 
Using xq = 4.4, and the result (34) for A, we find d = 
6 x 10~ 7 . We can compare (36) to that derived from the 
exact solution (27) of the Starobinsky equation. In this 
case we find 



d. 



hx\ 
yTj2 



12 



1/2 



T 2 

x . 



(37) 



Using the result (35) for A s , and x = 4.4, we find 
d s = 5 x 10~ 7 . We therefore see that over scales of the ob- 
servable universe both solutions predict nearly identical 
and negligable deviations from gaussian. 

By assuming xq = xqq we are simply trying to ex- 
clude the effects of fluctuations on scales larger than our 
present horizon. In actual fact inflation would most likely 
have been in progress for a long time before the observ- 
able universe left the Hubble radius during inflation. In 
this case we can consider x to be anywhere in the range 
4.4 < x < x p i. For the usual self-coupling constant in 
(35) we find that the Planck boundary is x pl ~ 5000. 
Because the deviations from gaussian in (37) are propor- 
tional to Xq there is a possibility of significant deviations 
from gaussian fluctuations on super-horizon scales. For 
the new coupling constant in (34) we have x p i ~ 24. 
In this case the deviations from gaussian in (36) depend 
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only linearly on xq . This combined with the much smaller 
Planck boundary means that the solution (23) predicts 
gaussian fluctuations to a high accuracy over both super 
and sub-horizon scales. We have ignored the effects of 
boundary conditions that should be imposed at <p e and 
4> p i. These effects will be small and and are unlikely to 
change these conclusions. 

Yi, Vishniac and Minehige [7] have analysed in de- 
tail solutions of the Starobinsky equation (in both Ito 
and Statonovich interpretations) for the model discussed 
here. They discussed deviations from gaussian in terms 
of a simple measure for skewness derived from the prob- 
ability distribution of the exact solution (27). We have 
also applied this measure to the exact solution (23) and 
found it to give results consistent with the simpler mea- 
sure discussed above. 

In this letter we showed that the new theory of in- 
flaton dynamics developed in [16] can address the fine- 
tuning problem yet still predict density fluctuations in 
general agreement with observations. This was demon- 
strated for the quartic potential but is likely to gener- 
alize to a wide range of potentials. These results were 
based on the classical limit of the theory. The other 
great advantage of this theory is that it leads naturally 
to a description of the inflaton as a quantum open system. 
This allows the quantum-to-classical transition to occur 
as a non-equilibrium quantum statistical process (deco- 
herence), rather than being simply postulated as in the 
conventional approach. Work on this issue is currently 
in progress. 
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